Jets, Bulk Matter, and their Interaction in Heavy Ion Collisions at Several TeV 
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We discuss a theoretical scheme that accounts for bulk matter, jets, and the interaction between 
the two. The aim is a complete description of particle production at all transverse momentum (pt) 
scales. In this picture, the hard initial scatterings result in mainly longitudinal flux tubes, with 
transversely moving pieces carrying the pt of the partons from hard scatterings. These flux tubes 
constitute eventually both bulk matter (which thermalizes and flows) and jets. We introduce a 
criterion based on parton energy loss to decide whether a given string segment contributes to the 
bulk or leaves the matter to end up as a jet of hadrons. Essentially low pt segments from inside the 
volume will constitute the bulk, high pt segments (or segments very close to the surface) contribute to 
the jets. The latter ones appear after the usual flux tube breaking via q-qbar production (Schwinger 
mechanism). Interesting is the transition region: Intermediate pt segments produced inside the 
matter close to the surface but having enough energy to escape, are supposed to pick up q-qbar 
pairs from the thermal matter rather than creating them via the Schwinger mechanism. This 
represents a communication between jets and the flowing bulk matter (fluid-jet interaction). Also 
very important is the interaction between jet hadrons and the soft hadrons from the fluid freeze-out. 
We employ the new picture to investigate Pb-Pb collisions at 2.76 TeV. We discuss the centrality 
and Pt dependence of particle production and long range dihadron correlations at small and large 
Pt- 
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I. INTRODUCTION 

Traditionally the physics of ultrarelativistc heavy ion 
collisions is discussed in terms of different categories like 
collective dynamics, parton-jet physics, and fluctuation- 
correlation studies, although these different topics are 
highly correlated. In this article, a complete dynamical 
picture of particle production at all pt scales will be pre- 
sented, which accounts for the production and evolution 
of bulk matter and jets, and the very important inter- 
action between the two components (which is not only 
the well known parton energy loss). The consequences 
of these interactions can be nicely seen in long range di- 
hadron correlations. 

The physical picture of our approach is the following: 
Initial hard scatterings result in mainly longitudinal flux 
tubes, with transversely moving pieces carrying the pt of 
the partons from hard scatterings. These flux tubes con- 
stitute eventually both bulk matter (which thermalizes, 
flows, and finally hadronizes) and jets, according to some 
criteria based on partonic energy loss. We will consider a 
sharp fluid freeze-out hypersurface, deflned by a constant 
temperature. Freeze-out here simply means the end of 
the fluid phase, but the hadrons still interact. High en- 
ergy flux tube segments will leave the fluid, providing jet 
hadrons via the usual Schwinger mechanism of flux-tube 
breaking caused by quark-antiquark production. 

But the jets may also be produced at the freeze-out sur- 
face. Here we assume that the quark-antiquark needed for 



the flux tube breaking is provided by the fluid, with prop- 
erties (momentum, flavor) determined by the fluid rather 
than the Schwinger mechanism. Considering transverse 
fluid velocities up to 0.7c, and thermal parton momentum 
distributions, one may get a "push" of a couple of GeV 
to be added to the transverse momentum of the string 
segment. This will be a crucial effect for intermediate pt 
jet hadrons. 

There is another important issue. Even for hadrons 
with transverse momenta of 10-20 GeV, there is a large 
probability of a jet hadron formation before it enters the 
dense hadronic medium. This means a signiflcant proba- 
bility of scatterings between jet hadrons and soft hadrons 
(from freeze-out), having essentially two consequences: 
an increase of low pt particle production, and a reduction 
of yields at high pt . In addition there are of course the well 
known hadronic interactions between the soft hadrons. 

We have discussed different processes which all affect 
Pt spectra. It is, however, possible to disentangle dif- 
ferent contributions, by looking at dihadron correlations. 
These are extremely useful tools heavily used by exper- 
imental groups at the RHIC and the LHC [iHal- Re- 
cently, the CMS and ALICE collaborations published 
results on such correlations in Pb-Pb collisions at 2.76 
GeV and different centralities, over a more or less broad 
range in relative pseudorapidity (Att) and full coverage of 
the relative azimuthal angle (A;/)) [J, l5|. Different com- 
binations of transverse momenta pf^^°^ and p"^^ of as- 
sociated and trigger particles in the range between 0.25 
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GeV/c and 15 GeV/c are investigated. Considering long 
range correlations (IA77I > A, A > 0.8), the coeffi- 
cients VnA of the harmonic decomposition factorize as 
VnA = v{pf^^°'^) v{pf"^^) - not only for small transverse 
momenta but also for example for large p"^^ and small 
passoc_ Pqj, guY^w momenta the situation seems to be clear: 
the correlation is flow dominated. But factorization does 
not necessarily mean that both hadrons carry the flow 
from the fluid! This can in particular not be the answer 
for observed correlations for large p^^^^ - here we have 
to deal with an interaction between the flowing bulk and 
jets, which makes the observed correlations very interest- 
ing in particular as a test of our ideas concerning bulk-jet 
separation and interaction. 

Another challenge : The ATLAS collaboration showed 
recently results [6| on elliptical flow of charged particles 
with respect to an event plane in the opposite 77 hemi- 
sphere (also a kind of long range correlation). The V2 
values are quite large up to values of pt — 20GeV/c, for 
eight different centrality ranges. Can we understand this 
in a quantitative fashion? 

The heavy ion results shown in this paper are based 
on 2000000 events simulated with EPOS2.17v3. A cen- 
tral (0-5%) Pb-Pb event takes on the average around 2 
HS06 hours CPU time, using machines with an average 
scaling factor of 8.7 Q. Always six events share the same 
parton configuration and hydrodynamic evolution, with 
only particle production and hadronic rescattering being 
redone (to gain statistics). This is taken care of when 
considering mixed events in correlation studies. 



II. THE BASIS: FLUX TUBES FROM A 
MULTIPLE SCATTERING APPROACH 



The starting point is a multiple scattering approach 
corresponding to a marriage of Gribov-Regge theory and 
perturbative QCD (pQCD), see Fig. [TJ An elementary 
scattering corresponds to a parton ladder, containing a 
hard scattering calculable based on pQCD, including ini- 
tial and final state radiation (for details see [8]). These 
ladders are identified with flux tubes (see Fig. [5]), which 
are mainly longitudinal objects, with transversely mov- 
ing parts, carrying the transverse momenta of the hard 
scatterings. These objects are also referred to as kinky 
strings. One should note that here multiple scattering 
does not mean just a rescattering of hard partons, it 
rather means a multiple exchange of complete parton lad- 
ders, leading to many flux tubes. In this case, the energy 
sharing between the different scatterings will be very im- 
portant, to be discussed later. 

The consistent quantum mechanical treatment of the 
multiple scattering is quite involved, it is based on cutting 
rule techniques to obtain partial cross sections, which are 
then simulated with the help of Markov chain techniques 




Figure 1: (Color online) Multiple scattering in nucleus- 
nucleus collisions. 




Figure 2: (Color online) An elementary parton ladder, whose 
final state is identified with a color fiux tube (kinky string). 



i- 

As said before, the final state partonic system cor- 
responding to elementary parton ladders are identified 
with flux tubes. The relativistic string picture [l(]|4l2| 
is very attractive, because its dynamics is essentially de- 
rived from general principles as covariance and gauge in- 
variance. The simplest possible string is a surface A (a, /3) 
in 3+1 dimensional space-time, with piecewise constant 
initial velocities dX/df3. These velocities are identified 
with parton velocities, which provides a one to one map- 
ping from partons to strings. For details see [3>@|- The 
high transverse momentum (pt) partons will show up as 
transversely moving string pieces, see Fig. ^a.). Despite 
the fact that in the TeV energy range most processes are 
hard, and despite the theoretical importance of very high 
Pt partons, it should not be forgotten that the latter pro- 
cesses are rare, most kinks carry only few GeV of trans- 
verse momentum, and the energy is nevertheless essen- 
tially longitudinal. In case of elementary reactions, the 
strings will break (see Fig. ^h) via the production of 
quark-antiquark pairs according to the so-called area law 
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Figure 3: (Color online) (a) Flux tube with transversely mov- 
ing part (kinky string) in space, at given proper time, 
(b) Flux tube breaking via q — q production, which screens 
the color field (Schwinger mechanism). 



[a, 0, [ij, [Tj| . The string segments are identified with final 
hadrons and resonances. 

This picture has been very successful to describe par- 
ticle production in electron-positron annihilation or in 
proton-proton scattering at very high energies. In the 
latter case, not only low pt particles are described cor- 
rectly, for example for pp scattering at 7 TeV [isl [l6fl . 
also jet production is covered. As discussed earlier, the 
high transverse momenta of the hard partons show up 
as kinks, transversely moving string regions. After string 
breaking, the string pieces from these transversely mov- 
ing areas represent the jets of particles associated with 
the hard partons. To demonstrate that this picture also 
works quantitatively, we compute the inclusive pt distri- 
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Figure 4: (Color online) Inclusive pt distribution of jets. We 
show the calculation (full lines) compared to ATLAS data ,18] 
(dotted lines). 
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Figure 5: (Color online) Transverse momentum distribution 
of partons: We compare our results (red stars) with a parton 
model calculation (dotted line). 



bution of jets, reconstructed with the anti-kt algorithm 
(l7J and compare with data [Hi, see Fig. g) In Fig. O 
we compare our pt distribution of partons with a parton 
model calculation based on CTEQ6 parton distribution 
functions [1^, in both cases leading order with a K- factor 
of 2. 



III. JET BULK SEPARATION 

In heavy ion collisions and also in high multiplicity 
events in proton-proton scattering at very high energies, 
the density of strings will be so high that the strings can- 
not decay independently as described above. Here we 
have to modify the procedure as discussed in the follow- 
ing. The starting point are still the flux tubes (kinky 
strings) originating from elementary collisions. These 
flux tubes will constitute both, bulk matter which ther- 
malizes and expands collectively, and jets. The criterion 
which decides whether a string piece ends up as bulk or 
jet, will be based on energy loss. In the following we con- 
sider a flux tube in matter, where "matter" first means the 
presence of a high density of other flux tubes, which then 
thermalize. A more quantitative discussion will follow. 

Three possibilities should occur, referred to as A, B, C, 
see Fig. [H] : 

A String segments far from the surface and/or be- 
ing slow will simply constitute matter, they loose 
their character as individual strings. This matter 
will evolve hydrodynamically and finally hadronize 
("soft hadrons"). 

B Some string pieces (like those close to transversely 
moving kinks) will be formed outside the matter, 
they will escape and constitute jets ("jet hadrons"). 

C There are finally also string pieces produced inside 
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matter or at the surface, but having enough en- 
ergy to escape and show up as jets ("jet hadrons"). 
They are affected by the flowing matter ("fluid-jet 
interaction") . 

Let us discuss how the above ideas are reahzed. In prin- 
ciple the formation and expansion of matter and the in- 
teraction of partons with matter is a dynamical process. 
However, the initial distribution of energy density and the 
knowledge of the initial momenta of partons (or string 
segments) allows already an estimate about the fate of 
the string segments. By "initial time" we mean some early 
proper time tq which is a parameter of the model. Strictly 
speaking, energy loss concerns partons, modifying even- 
tually the kink momenta in our picture, and the momenta 
of the string segments after breaking will be reduced. We 
will therefore base our discussion on energy loss on string 
segments. 

We estimate the energy loss AE of string segments 
along their trajectory to be 



AE = fcEioss Eo / (pK))'/' max(l, ^E/Eq ) dL/Lo, (1) 



inspired by [20] , where p is the density of string segments 
at initial proper time tq , Vq is an elementary volume cell 
size (technical parameter, taken to be 0.147 fm^), Lq is 
a (technical) length scale (taken to be 1 fm), E the en- 
ergy of the segment in the "Bjorken frame" moving with 
a rapidity y equal to the space-time rapidity ris, dL is a 
length element, and fegioss and Eq are parameters. We 



(a) 






(b) 




Figure 6: (Color online) Flux tube in matter (from other flux 
tubes, blue colored area). One distinguishes three types of 
behavior for string segments, noted as A, B, C (see text). 
The highest pt string segment may be (a) of type B or (b) of 
type C. 



introduce an energy cutoff Eq to have sufHcient energy 
loss for slowly moving segments. 

A string segment will contribute to the bulk (type A 
segment), when its energy loss is bigger than its energy, 
i.e. 



AE> E. 



(2) 



All the other segments are allowed to leave the bulk (type 
B or C segments). Only the bulk segments are used 
to determine the initial conditions for hydrodynamics, 
following the same procedure as explained in [8| (with 
some new elements, as discussed in the next section). 
Starting from this initial condition, the bulk matter will 
evolve according to the equations of ideal hydrodynam- 
ics till "hadronization", which occurs at some "hadroniza- 
tion temperature" Th [8]. Hadronization means that we 
change from matter description to particle description, 
but hadrons still interact among each other, realized via 
a hadronic cascade procedure [21], already discussed in 
[8,]. _ 

After having performed the hydrodynamic expansion, 
we have to come back to the string segments which escape 
the bulk because their energy is bigger than the energy 
loss. We employ a formation time: the string segments 
are formed at times t distributed as exp(— i/7rform); with 
some parameter rform which is taken to be Ifm/c. If the 
formation time is such that the segment is produced out- 
side the "hadronization surface" defined by Th, the seg- 
ment will escape as it is (type B segment). 

Most interesting are the segments which are formed in- 
side but still escape, because they have E > AE. These 
are type C segments. They escape, but their properties 
change. Actually such a segment leaves "matter" at the 
hadronization surface at a particular space-time point x, 
which is characterized by some collective flow velocity 
v{x). We assume that the string breaking in this case is 
modified such that the quark and antiquark (or diquark) 
necessary for the string breaking are taken from the flow- 
ing fluid rather than being produced via the Schwinger 
mechanism. So the new string segment is composed of 
a quark and antiquark (diquark) carrying the flow veloc- 
ity, and the string piece in between, which has not been 
changed. This string piece may or may not carry large 
momentum, depending on whether it is close to a kink or 
not, the former possibility shown in Fig. [71 

In any case, due to the fluid-jet interaction, the prop- 
erties of this segment change drastically compared to the 
normal fragmentation: 

• The quark and antiquark (or diquark ) from the 
fluid provide a push in the direction of the moving 
fluid. 

• The quark (antiquark) flavors are determined from 
Bose-Einstein statistics, with more strangeness pro- 
duction compared to the Schwinger mechanism. 
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PbPb 2.76 TeV 



Figure 7: (Color online) A type C segment picks up quark 
and antiquark from the fluid, carrying momenta and flavor 
according to the fluid properties ("fluid-jet interaction"). 



• The probability p^iq to have a diquark rather 
than an antiquark will be bigger compared to a 
highly suppressed diquark-antidiquark breakup in 
the Schwinger picture (pdiq is a parameter). 

Our procedure has 4 parameters: /ceioss (=0.042), Eq 
(=6 GeV), Tform(=l fm/c), pdiq(=0.22). It allows to cover 
in a single scheme the production of jets, of bulk, and the 
interaction between the two. 



IV. FORMATION TIMES 

A crucial ingredient to the mechanism of fluid-jet inter- 
action is the formation time of jet hadrons (the hadrons 
which leave the fluid) . The probability distribution of the 
formation times t of jet hadrons with gamma factors 7 is 
given as 



prob ^ exp 



t 



(3) 



7'''form , 

(we use Tform — Ifm/c). The probability of having a 
formation point inside the fluid is obtained as an integral 
over eq. ([3]), 

1-expf-^^^), (4) 

with tniax being the time corresponding to a formation at 
the fluid surface. Rather than making a simulation, we 
are going to present a very simple formula providing a 
rough estimate of the pt dependence of this probability. 
For a collision of two Pb nuclei in some centrality interval, 
characterized by the mean impact parameter 6, we use 
ctmax = fph — b/2, whcrc rpb is the radius parameter used 
in the Wood-Saxon distribution of nucleons. Considering 
transversely moving hadrons of mass m, we have 7 w 
Pt/mc. The estimate Pinsidc of the probability to form 
(pre)hadrons inside the fluid is 




exp 



(''Pb - b/2) m 



PtTio 



(5) 



" 5 10 15 20 25 30 35 40 45 50 
P, (GeV/c) 

Figure 8: (Color online) The estimate Pinsido to form 
(pre)hadrons inside the fluid, as a function of pt, for Pb-Pb 
collisions at 2.76 TeV. We show the curves for the 0-5% and 
the 20-30% most central events. 



In Fig. [H we show the result for the 0-5% and the 20- 
30% most central events in Pb-Pb collisions at 2.76 TeV, 
using CTform = 1 fm, mc^ = 1 GeV, rpb = 6.5 fm, and 
for the average impact parameters b = 1.8 fm (0-5%) and 
& = 7.8fm (20-30%). 

By construction, the probability Pfluid-jct of having a 
fluid-jet interaction is equal to the probability of forming 
(pre)hadrons inside the fluid, so its estimate is given by 
^inside- From Fig. [51 we see that the probability is quite 
large for intermediate values oi pt, but even large values 
(50 GeV/c) are significantly affected. Whether the effect 
of the interaction can be seen in some observable is a 
different question and will be discussed later. 

Several authors have already discussed about "in- 
medium hadronization", see for example Ref. |22| , where 
one also finds an overview about earlier models on this 
subject. 



V. HYDRODYNAMICS 

The bulk matter extracted as described above provides 
the initial condition for a hydrodynamic evolution. As ex- 
plained in [Sj, we compute the energy momentum tensor 
and the flavor flow vector at some position x (at r = tq) 
from the four-momenta of the bulk string segments. The 
time T = To is as well taken to be the initial time for the 
hydrodynamic evolution. This seems to be a drastic sim- 
pliflcation, the justification being as follows: we imagine 
to have a purely longitudinal scenario (descibed by flux 
tubes) till some proper time thux < '''o • During this stage 
there is practically no transverse expansion, and the en- 
ergy per unit of space-time rapidity does not change. This 
property should not change drastically beyond taux, so we 
assume it will continue to hold during thermalization be- 
tween Tfiux and tq. So although we cannot say anything 
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Figure 9: (Color online) Energy density and pressure versus 
temperature, for our equation-of-state (lines) compared to lat- 
tice data [23| (points). 



about the precise mechanisni which leads to thermaliza- 
tion, and therefore we cannot compute the real T^'', we 
expect at least the elements T^'^ and T°' to stay close to 
the flux tube values, and we can use the flux tube results 
to compute the energy density. Only T*^ will change con- 
siderably, but this does not affect our calculation much. 
We employ three-dimensional ideal hydrodynamics as 
described in Q, with some modification to be discussed 
in the following. As in [8|, we construct the equation of 
state as 



p = PQ + X{pH -Pq), 



(6) 



where pn is the pressure of a resonance gas, and pg the 
pressure of an ideal quark gluon plasma, including bag 
pressure. We use an updated A: 



A = exp {-z - 3z') e(T - T,) + Q{T, - T), (7) 



with 



z = .t/(1 + .t/0.77). 



{T-T,)/6, 



(8) 



using S = 0.24exp(— /i^/0.4^). The new A provides an 
eqtiation of state in agreement with recent lattice data 
dil, see Fig. [i 

Apart of the new equation of state, we use the same 
procedure to obtain energy density and pressure from the 
string segments, as described in Q. However, 

• doing the calculation for Pb-Pb collisions at 2.76 
TeV, we get too much elliptical flow (20-30%), a 
hint that one should include viscosity. 

Taking the usual small radii of the elementary flux-tubes, 
we get extremely strongly fluctuating energy densities (in 
the transverse plane) . Viscosity will quickly reduces these 
strong fluctuations. 



200 



150 



100 




Figure 10: (Color online) Initial energy density in a central 
Pb-Pb collision at 2.76 TeV, at a space-time rapidity rjs = 0, 
as a function of the transverse coordinates x and y. We take 
artificially large values of the flux tube radii - which provides 
relatively smooth initial conditions - to mimic viscous effects. 



• We try to mimic viscous effects by taking artificially 
large values of the flux tube radii (we take 1 fm), 
in order to get smoother initial conditions. 

This has the effect of reducing the elliptical flow by 20- 
30%, as needed. In Fig. [TUl we show an example of such 
an initial energy density. 



VI. CENTRALITY DEPENDENCE OF THE 
MULTIPLICITY AND FAKE SCALING LAWS 

As a very first check of the new approach, we consider 
the centrality dependence of the charged particle yield. 
Although very basic, there is quite some confusion about 
this quantity. Whereas hard processes scale roughly with 
the number of binary nucleon-nucleon collisions (in a sim- 
ple geometrical picture), the centrality dependence of the 
charged particle yield (dominated by low pt) is very dif- 
ferent: it looks more like a scaling with the number of 
participating nucleons. This reminds us of the good old 
"wounded nucleon model", which has a physical mean- 
ing - at low energies: the projectile and target nucleons 
are excited ("wounded"), and this is the main source of 
particle production. 

Amazingly, this approximate participant scaling holds 
also at higher energies, the centrality dependence at the 
LHC is almost identical to the one at the RHIC [24] . This 
is quite strange, since one might believe that at higher 
energies hard processes dominate, so one could expect 
more binary scaling. But this is not the case. 
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What do we 
proach? In Fig. 



get in the multiple scattering ap- 
1111 we plot the yield per participant 




10 20 30 40 50 60 70 80 90 100 

centrality (%) 

Figure 11: (Color online) Centrality dependence of the 
charged particle yield. Calculations (stars) are compared to 
ATLAS data [H (circles). 

( dn/dr]{Q)/Npart ) as a function of the centrality. Npart 
is the number of participating nucleons. As in the data 
[25j , we obtain a moderately increasing yield per partici- 
pant. How can this happen? How can one get something 
like a wounded nucleon result at the LHC? In the model 
we can of course easily check the relative contribution of 
particle production from remnant decays. In Fig. 1121 we 
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Figure 12: (Color online) Relative fraction of particle produc- 
tion from remnants versus rapidity. The contribution from 
remnants at mid-rapidity is very small. 

plot the relative fraction of particle production from rem- 
nants as a function of the rapidity. As expected, remnant 
particle production is important at large rapidities, but 
the contribution at mid-rapidity is close to zero. So the 
physical mechanism of soft particle production is not a 
wounded nucleon picture. 

In our approach, the source of particle production is 
the flux tubes, originating from elementary scatterings, 
which are in principle proportional to the number of bi- 
nary nucleon-nucleon collisions. But there are impor- 




Figure 13: (Color online) Energy sharing in AA scatterings. 
The sum of all "ladder end" light cone momentum fractions 
a; J, linked to a given remnant and the remnant fraction x^^^^^ 
have to add up to unity. 



tant effect due to energy conservation and shadowing, 
discussed in detail in [^ |^. In our multiple scattering 
approach (which determines the initial conditions), the 
complete AA scattering amplitude is expressed in terms 
of elementary contributions, which are parton ladders, 
later showing up as strings. Each parton ladder is char- 
acterized by the light cone momentum fractions x\ and 
x^ of the "ladder ends", which are the outer partons of 
the ladder, see Fig. [T3](also transverse momenta are con- 
sidered, but not discussed here). 

It is a unique feature of our approach that we do a pre- 
cise bookkeeping of energy and momentum: For each nu- 
cleon (projectile or target) the initial energy-momentum 
has to be shared by all the ladders connected to this nu- 
cleon and the nucleon remnants, i.e. for all nucleons i we 
have 



E 



■-fc 



1, (9) 



all ladders k 
connected to nucleon i 



where x^^^^^ is the momentum fraction of the nucleon 
remnant i. These are very strong conditions, which affect 
the results substantially, see [9|. 

The most important consequence relevant for our dis- 
cussion here is the fact that parton ladders leading to low 
Pt particles are suppressed compared to what is expected 
from binary scaling. We will get a nuclear modification 
factor which is less than one at low pt ■, as shown in Fig. [TJ] 
for central Pb-Pb collisions at 2.76 TeV. So although par- 
ticle production at central rapidities in very high energy 
collisions is dominated by binary scattering (providing 
the initial energy density), particle production does not 
increase proportional to the number of binary scatterings, 
due to energy conservation. 

It is absolutely necessary that binary scaling is broken 
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Figure 14: (Color online) Nuclear modification factor Raa = 
N~^l^{dnAA/dpt) / (dupp/dpt) vs transverse momentum pt, 
showing the breaking of binary scaling at low pt (due to en- 
ergy conservation). The resulting "approximate participant 
scaling at low pt" is a pure coincidence! 



at low Pt, because it is simply an experimental fact. The 
usual explanation is a two component picture: hard scat- 
tering at high Pt which shows binary scaling and a soft 
component which scales as the number of participants In 
our picture, binary collisions determine everything. But 
certain binary coUisions are suppressed due to energy con- 
servation, leading to a deviation from Raa — 1- 

We will discuss the pt dependence of Raa in the next 
section. Here we present for completeness the pseudora- 
pidity distributions of charged particles for different cen- 
tralities, see Fig. 1151 where we compare our calculation 
with data from ATLAS Ei. 



Figure 15: (Color online) Pseudorapidity distributions of 
charged particles for different centralities. The lines are cal- 
culations, the the circles are ATLAS data [25|, see text. 
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VII. TRANSVERSE MOMENTUM 

DEPENDENCE OF PARTICLE YIELDS: 

IMPORTANCE OF HADRONIC RESCATTERING 

OF SOFT AND JET HADRONS 

We first investigate particle production at low trans- 
verse momenta. In figs. [1^1 [iTl and [THl we show trans- 
verse momentum distributions of pions, kaons, and pro- 
tons, for central and semi-peripheral Pb-Pb collisions at 
2.76 TeV. We compare the full calculation including hy- 
drodynamic evolution and hadronic final state cascade 
(solid lines) with the calculation without cascade (dashed 
lines) and with data from ALICE {20|. 

In order to understand the results, one has to recall 
that not only the "soft" particles produced from the fluid 
may interact, but also the jet-particles having enough en- 
ergy to escape the fluid may interact with these soft par- 
ticles. In particular intermediate pt jet particles are can- 



Figure 16: (Color online) Transverse momentum distributions 
of (from top to bottom) negative pion, kaons, and protons, in 
the 0-5% most central Pb-Pb collisions at 2.76 TeV. We show 
the full calculation (solid lines) and the ones without hadronic 
cascade (dashed lines), compared to ALICE data (circles) |2(j |. 



didates, because their formation time will produce them 
just in the high density hadronic region. Let us discuss 
the consequences of these interactions, by comparing the 
solid and dashed curves in the figures. 

• We see in particular in Fig. \W\ a strong reduction 
of protons at low pt due to hadronic rescattering, 
which can be attributed to proton-antiproton anni- 
hilation among the soft hadrons. 

• We see also a sizable increase of pion production at 
low Pt, which is due to inelastic rescatterings of jet 
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Figure 17: (Color online) Same as Fig. \T^ but 20-30% most 
central. 
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Figure 18: (Color online) Same as Fig. [161 but 40-50% most 
central. 



hadrons with soft ones 

In figs, [ini [T71 and [THl we only show results up to 3 
GeV/c, because this is the range where data on protons, 
pions, and kaons are available. It is nevertheless inter- 
esting to know the effect of jet-soft scattering beyond 3 
GeV/c. We therefore plot in Fig. ^^ the ratio of the 
full calculation to the one without hadronic cascade, for 
the pt spectra of charged particles (dominated by pions) 
in central Pb-Pb coUisions at 2.76 TeV, up to 20 GeV/c. 
There is a big effect at intermediate values of pt - up to 
20 GeV/c ! In other words, jet-soft rescattering is very 
important in this range. Similar observation have already 
been made in [23| for AuAu collisions at the RHIC. 
The big effect of the jet-soft interaction can be under- 



Figure 19: (Color online) Transverse momentum dependence 
of the ratio of the full calculation to the one without hadronic 
cascade, for charged particle production in the 0-5% most 
central Pb-Pb collisions at 2.76 TeV. 
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Figure 20: (Color online) Transverse momentum dependence 
of 1 — 7? (red dots), with R being the ratio plotted in Fig. 
1191 The dashed line is the estimate Pinsidcof the probability 
to produce a jet hadron inside the fluid. 



stood by plotting 1 — R, with R being the ratio (with / 
without cascade) plotted in Fig. [121 together with the 
probability estimate Pinsido to produce a jet (pre)hadron 
inside the fluid, see fig. [501 These early produced hadrons 
go through the dense hadronic phase (of soft hadrons), 
and Pinsidc is therefore also a measure of the probability 
of having a jet-soft interaction. We see indeed 



l^R^Pir 



(10) 



at large pt (in absolute terms, without adding factors). 
Even though we are running out of statistics, it is clear 
from the above discussion that the effect goes well beyond 
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Figure 21: (Color online) The nuclear modification factor in 
Pb-Pb at 2.76 TeV vs pt: We compare data [28] (squares) with 
the full calculation (red line + circles) and its jet contribution 
(yellow line + triangles) , as well as the bulk (hydro) contri- 
bution of a calculation without hadronic cascade (blue line + 
squares) . 



20 GeV/c. 

To compare the pt spectra with experimental data, one 
uses often the so-called nuclear modification factor Raa, 
which is the ratio of the inclusive transverse momentum 
spectrum of particles in nucleus-nucleus scatterings over 
the proton-proton ones, normalized by the number A^coU 
of binary collisions. Doing this procedure, we obtain the 
curves shown in Fig. I21[ where we plot our simulation 
results for charged particle production together with the 
data from ALICE [28]. We show the full model, includ- 
ing hydrodynamic evolution and final state hadronic cas- 
cade [21j and its jet contribution from the string segments 
which escaped from the bulk and which did not rescat- 
ter. We also show the bulk contribution (originating from 
the hadronized fiuid) from the calculation without final 
state hadronic cascade. The two latter curves do not add 
up to give the full result - the difference is due to the 
"secondary interactions" discussed earlier: 

• Fluid-jet interaction - pushing the jet hadrons at 
intermediate pt to higher transverse momenta. 

• Jet-soft interactions between jet hadrons and soft 
ones from fiuid freeze-out. 

There are also soft-soft interactions (among soft hadrons 
from fiuid freeze-out), which are important for baryon 
yields, but not so much for the charged particle results. 
From the above discussion it is clear that even con- 
sidering elementary quantities as charged particle yields, 
it is difficult to make any quantitative analysis without 



jet-soft -> yield increase 

soft-soft -> baryon annihilation 

soft -soft -> softer spectrum 

j et-soft -> softer s pectrum 

parton en. loss -> softer spectrum 

fluid-jet push -> harder spectrum 
-^ ^- 

1 2 To 20 pt(GeV/c) 



Figure 22: (Color online) Secondary interactions. The red 
arrows indicate the pt range which are affected. 



considering these "secondary interactions". We sketch the 
different interactions in Fig. [211 



VIII. DIHADRON CORRELATIONS IN PB-PB 
AT 2.76 TEV 

Our prescription for bulk-jet separation and interac- 
tion should also strongly affect dihadron correlations, 
which provide much more information than simple spec- 
tra. With all parameters (/ceioss, -E-q, Tseg, Pdiq) being 
fixed from the considerations in the last section, we now 
compute dihadron correlation functions defined as 

where S is the number of pairs in real events, and M the 
number of pairs for mixed events. 

As an example, we show in Fig. [53] the correlation 
function for the pt of the trigger particle {p"^^) in the 
interval 4.5-5.5 GeV/c and the pt of the associated parti- 
cle (passoc-) jj^ ^i^g j.j^j^gg 2_2.5 GeV/c, in the 0-10% most 

central Pb-Pb collisions at 2.76 TeV. Besides the jet peak 
at Acj) = and Arj = 0, we clearly identify a completely 
fiat ridge over the full range in Arj at A<j) — 0. 

The reason for the ridge structure is the fact that there 
is an azimuthal asymmetry of the initial energy density 
(see Fig. [TU)l . Although the energy density is biggest 
around space-time rapidity rjs = and drops fast towards 
forward and backward rjs , the shape of the asymmetry is 
preserved. This leads finally to an asymmetric flow, again 
very similar at different values of r]s , and this "makes" the 
long range correlation at A<j) — 0. 

The smooth rjs dependence of the energy density in our 
approach (see Fig. IM)) is due to the fact the energy den- 
sity is calculated from flux tubes. And these flux tubes 
have to be treated correctly as continuous longitudinal 
objects (as we do). In an earlier version, we treated flux 
tubes via the randomly (in r]s) distributed fiux tube seg- 
ments, obtained from a string fragmentation procedures. 
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This gives a bumpy structure in 77^- the ridge is not flat 
any more but has a Gaussian shape! So the flux-tube ba- 
sis is an essential ingredient for obtaining a perfect ridge 
shape, as observed in the data. 

In Fig. [551 we show a correlation function for p"^^ in 
the interval 5.5-8.0 GeV/c and pf"""" in the range 2-2.5 
GeV/c, again in the 0-10% most central Pb-Pb collisions 
at 2.76 TeV. Although the trigger pt is too large to orig- 
inate from freeze-out (from the flowing fluid), one still 
observes a ridge structure, which is due to the fluid-jet 
interaction. Let us consider again the situation of an ini- 
tial azimuthal anisotropy in the energy density which is 
transported into a corresponding anisotropy in the flow, 
as discussed earlier. We sketch in Fig. [Ulthe (somewhat 
exaggerated) situation of a triangular transverse flow pat- 
tern with maximal flow around = 0°, 120°, and 240° 
(with respect to the y-axis). The flow maxima are indi- 
cated by blue arrows. Again it is very important that 
this flow pattern is (not necessarily in magnitude, but in 
shape) very similar at different longitudinal positions - in 
the figure indicated by the two transverse planes P and 
P', corresponding to two different space-time rapidities 
rjs and ry^. A soft hadron (S) produced at rjs at the fiuid 
surface close to the position of maximal fiow (for exam- 
ple at (/) = 0°), will be boosted by the latter one and 
therefore carry information about this flow. A jet hadron 
(J) produced at 77^ at the same angle ((/> = 0°) close to 
the surface, will pick up a quark and an antiquark, both 
carrying flow, which adds the corresponding transverse 
momentum to the pt of the string segment (red element 
in the figure). It is the same flow which affects the jet 
hadron at r/^. and the soft hadron at rjs , which creates the 
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Figure 23: (Color online) Calculated dihadron correlation 
function for pf '^^ in the interval 4.5-5.5 GeV/c and pj^^^o'^ 
in the range 2-2.5 GeV/c in the 0-10% most central Pb-Pb 
collisions at 2.76 TeV. 
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Figure 24: (Color online) The energy density for a single event 
in a central Pb-Pb collisions at 2.76 TeV as a function of 
the longitudinal variable 77s and the transverse one y (and 
for a; = 0). 
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Figure 25: (Color online) Calculated dihadron correlation 
function for pf'*^«^ in the interval 5.5-8.0 GeV/c and pf''°'= 
in the range 2-2.5 GeV/c in the 0-10% most central Pb-Pb 
collisions at 2.76 TeV . 



dihadron correlation at Aip = 0, the "ridge". The correla- 
tion remains visible, even when the flow contribution to 
the jet hadron is only 10%, this is why the correlation is 
still present even for trigger transverse momenta beyond 
10 GeV/c. 

We will now discuss some examples of semi-peripheral 
Pb-Pb coUisions at 2.76 GeV/c. In Fig. ^ we show 
the correlation function for Pj"^^ in the interval 3.5- 
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Figure 26: (Color online) Sketch of two cuts of the fluid vol- 
ume corresponding to the space-time rapidities rjs and r;!, , the 
two corresponding transverse planes being P and P'. We show 
the example of a triangular flow pattern - the same at 773 and 
ri's- 
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Figure 28: (Color online) Calculated dihadron correlation 
function for pf'"^^ in the interval 5.5-8.0 GeV/c and pT""" 
in the range 2-2.5 GeV/c in the 40-50% most central Pb-Pb 
collisions at 2.76 TeV. The jet peak has been cut for better 
visibility. 



the jet peak, we clearly see a very similar elliptical flow 
structure as in the previous example. 

The correlation functions are essentially flat as a func- 
tion of Arj, for large Aij. One therefore gets complete 
information about the long range correlations by inte- 
grating over A77, 



i?(A0) 



1 



2(-B - ^) JA<\Ari\<B 



i?(Ar/, A(/)) dA?7, (12) 



where we use A — 0.8 and the maximum B = 2. This 
function agrees perfectly with its Fourier decomposition, 



Figure 27: (Color online) Calculated dihadron correlation 
function for pf'^^ in the interval 3.5-4.5 GeV/c and pf""" in 
the range 1-1.5 GeV/c in the 40-50% most central Pb-Pb col- 
lisions at 2.76 TeV. 



4.5 GeV/c and p^ssoc ^^ ^^^ range 1-1.5 GeV/c, in the 
40-50% most central Pb-Pb collisions at 2.76 TeV. It 
can be clearly seen from the figure that the elliptical 
flow {^ cos(2A0)) is dominant, besides the jet peak at 
A77 = 0, A(f) = 0. But also here higher order harmonics 
(^ cos(iA(/))) contribute, as we will discuss later. In Fig. 
1281 we show another example of a correlation function for 
semi-peripheral collisions, with somewhat bigger trigger 
pt. We use pt"sg in ^Yie interval 5.5-8.0 GeV/c and pf''°'= 
in the range 2-2.5 GeV/c. The Arj range is chosen smaller 
to avoid to big statistical fluctuations. The jet contribu- 
tion becomes dominant in this case, but when we cut off 



R{A(f)) = 1 + X! 2^"^ cos{nA(f>), 



(13) 



n=l 



using the first five terms. This is very convenient, because 
it allows to discuss the features of the correlation func- 
tions for different options for p"^^ and p'^^^"^ by simply 
considering the Fourier coefficients. 

In figs. [551 and [301 we plot some coefficients T/rtA as 
a function oi p"^^ for different intervals oi pf'^^°'^. The 
value of Pf"^^ is actually the mean value in a certain inter- 
val, the largest interval being 8-15 GeV/c. We compare 
our simulation (stars) with the results from ALICE Q 
(circles) . 

In the semi-peripheral collisions of Fig. [5^1 we see 
clearly the dominance of elliptical flow: the n — 2 co- 
efficients are by far the largest. Nevertheless, also the 
higher harmonics contribute. We see in all cases an in- 
crease of the coefficients with p'^^'^°^ and with p"^^ up to 
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Figure 29: (Color online) The Fourier coefficients Kia as 
a function of pf'^'^ for pT""" within 0.25-0.5 GeV/c (a), 1- 
1.5 GeV/c (b), 2-2.5 GeV/c (c), in the 40-50% most central 
Pb-Pb collisions at 2.76 TeV. We compare the ALICE data 
[3 (circles) with calculations (red stars) . 



values of around 2-3 GeV/c. At the latter values the hy- 
drodynamic flow contributes the most to the correlation 
between soft hadrons from the fluid. 

For higher transverse momenta, the coefficients get 
smaller, because the correlation between soft particles 



Figure 30: (Color online) Same as Fig. [29l but for the 0-10% 
most central Pb-Pb collisions at 2.76 TeV. We compare the 
ALICE data Q (circles) with calculations (red stars) . 



dies out. But V2Adoes not at all drop to zero at high 
Pt because here the correlations between soft and jet par- 
ticles come into play - the jet particles which suffered 
a push by the fluid, as discussed earlier (fluid-jet inter- 
action). The fluid transfers at maximum few GeV/c of 
transverse momentum to the jet, but this is easily visible 
in the correlation (even at 20 GeV/c). 
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Figure 31: (Color online) The Fourier coefficients V„a as a 
function of p*"^*^ for pf"" within 2-4 GeV/c in the 0-5% most 
central Pb-Pb collisions at 2.76 TeV. We compare the data 4] 
(squares) with calculations (red stars). 



The results for semi-peripheral collisions are very ro- 
bust and depend little on model parameters. The most 
important ingredient is the elliptical initial shape on the 
energy density, given by the nuclear geometry. The ef- 
fects depend of course on the flow velocity at freeze-out, 
but this is not a parameter but itself a robust result (with 
a maximum around 0.7c). Finally the results depend on 
the jet formation time, which should be around 1 fm/c 
(the value we actually took without really attempting a 
fine-tuning). 

The V2A coefficients for central Pb-Pb collisions (Fig. 
150)1 are first of all smaller as compared to the semi- 
peripheral ones (note the different scales in figs. [29] and 
[5n)l . simply because the dominant effect of large initial 
ellipticity from the geometry is absent. Here the elliptic- 
ity is purely random. Apart from this, we observe the 
same features as for the semi-peripheral collisions: in- 
crease with transverse momenta up to 2-3 GeV/c, then 
decrease. A big difference in central Pb-Pb collisions com- 
pared to semi-peripheral ones is the fact that the higher 
harmonics and in particular V3A contribute substantially, 
because here both elliptical and triangular initial shape 
are of random origin (and therefore comparable) , whereas 
for more peripheral collisions the geometrical elliptical 
shape dominates everything else. 

In Fig. [3D]it seems that our calculation underestimates 
V2A, in particular for the largest p^''''°'^ range (2-2.5. GeV). 
Fortunately, similar data exist from CMS [1| , for p^ssoc ^^ 
the range 2-4 GeV/c in the 0-5% most central Pb-Pb col- 
lisions. In Fig. [211 we plot the corresponding coefficients 
V2A and VsA- Here we slightly overpredict the data! 



IX. V2 AND FORMATION TIMES 

Whereas dihadron correlations provide the most com- 
plete information about particle production - in partic- 
ular concerning the role of the "flowing" fluid, one may 
get the essential information by considering the elliptical 
flow coefficient V2 of single particle production, which is 
defined as 



V2 = (cOs[2((?!)-0Rof)]), 



(14) 



where is the azimuth angle of a particle, and 0Rof some 
reference plane. In [6], for the particles in the forward 
(backward) 77 hemisphere, the reference plane is the event 
plane angle ^backward ^-^forward^^ obtained from counting 
all particles in the opposite hemisphere. The angles are 
obtained from 



ibackward/ff 



'orward 



1 

= -tan 



(sin2(/)) 
(cos 20) ' 



(15) 



where the average is done in the forward / backward rj 
hemisphere within 3.2 < |ry| < 4.8. The 112 coefficient is 
then computed as 



V2 



i forward/backward "i 



)]> 



(16) 



Resolution correction is taken care of by dividing this ex- 
pression hy R ^ v/(cOs[2A</>]), with A0 = ^backward _ 

0forward_ ^^ jj-^ j,^£ ^gj j^elating particles with the event 
plane of the opposite hemisphere, we have kind of a long 
range correlation, but less clean than using dihadron cor- 
relations with a A77 > A requirement. But as mentioned 
before, the essential features can be seen as well. In Fig. 
1321 we plot V2 as a function of the transverse momentum 
for different centralities in Pb-Pb collisions at 2.76 TeV. 
The magnitude of the elliptical flow coefficients increase 
at low pt to reach a maximum around 2-3.5 GeV/c and 
then drop slowly at large pt . 

The behavior at high pt is the most interesting aspect: 
even at 10 GeV/c, there is a significant amount of ellip- 
tical flow, due to the fluid-jet interaction, which pushes 
jet particles in the direction of the collective flow at the 
freeze-out surface (and this effect will continue up to even 
higher pt, but we are simply running out of statistics). 
The high pt behavior is closely related to the formation 
time discussion we had earlier. The non- vanishing V2 at 
high Pt is mainly due to fluid-jet interactions, so the val- 
ues should be related to the estimated probability Pinside 
to form the jet hadron inside the fluid, which is equiva- 
lent to the fluid-jet interaction probability. In Fig. [331 we 
show Pinside (multiplied by an arbitrary factor) , together 
with the calculated and experimental V2 already show in 
Fig. [221 which shows that V2 is indeed proportional to 
the fluid-jet interaction probability. To compute Pinside 
according to eq. ([S]), we use CTform — 1 fm, itic^ — 1 GeV, 
rpb = 6.5 fm, and b = 11.5 fm (50-60%). 
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Figure 32: (Color online) pt dependence of elliptical flow (de- 
fined with respect to the opposite hemisphere sub-event plane) 
for different centralities in Pb-Pb collisions at 2.76 TeV. We 
compare the ATLAS data [y| (circles) with calculations (red 
lines) . 
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Figure 33: (Color online) Same data and model calculation as 
shown in Fig. |32]for the case of 50-60% most central collisions, 
together with Pinside. 



Assuming a flat (j) distribution, the probability distri- 
bution for ip is {2tt)^ ^ d(f> / dif) . which is in the case of 



soft 

1 + lLdCosW. 

Pt 



(19) 



which should only be considered for — 7r/2 < tp < tt/2, 
since ip and (?!)fiow have to correspond to the same hemi- 
sphere. We get anisotropics of the order of pf °^* /pj"'"'^ , 
which means at transverse momenta around 10 GeV/c, a 
soft push as little as 1 GeV/c can produce anisotropics of 
the order of 10%. 



Having a finite Pinside is a necessary condition to get a 
substantial V2 at some large value of pt , but not a suffi- 
cient one. Let us therefore estimate the effect of the fluid- 
jet interaction on V2 ■ One may consider a toy model, with 
soft particle emission due to flow into some preferred az- 
imuthal direction, say 0flow • Let us assume a jet hadron 
getting pushed by the fluid in the direction of (/)flow , which 
corresponds to adding some p/*"^* to the "hard" transverse 
momentum, p^^'^'^ of the flux-tube segment. Without loss 
of generality, we may set 0fiow = 0. The total transverse 
momentum of the jet hadron is 



;*Soft 



Pt 



■Pt 



d_ /^pf°"-|-pf^^^cos(0) 
p^''"^ sin(</>) 



ict 



cosV' 
simp 



(17) 
where ■0 is the jet hadron direction with respect to the 
flow direction (pfiow — 0. We have 



tan?/; = 



p^"^'^ sin((/>) 



psoh _!_ phard cos{(j)) ' 



(18) 



X. COALESCENCE 

For many years, different models have been employed 
to treat particle production at different transverse mo- 
mentum scales. The so-called intermediate range from 
2 to 6 GeV/c has been the domain of coalescence mod- 
els |29l - [33| . where hadrons are produced by recombining 
quarks from the plasma, to be distinguished from "frag- 
mentation" of partons. 

In our picture, there are certain aspects which give sim- 
ilar results as coalescence, but it is not a coalescence ap- 
proach. Already the notion "intermediate pi" extends to 
say 20 GeV/c and not 6. The corresponding transverse 
momentum of hadrons does not originate from plasma 
quarks and antiquarks - the main part is coming from 
the original flux tube. Whereas usual flux tube breaking 
in vacuum creates quark-antiquark pairs via a tunnel- 
ing process, the fluid-jet interaction amounts to replacing 
these quark-antiquark pairs by partons from the plasma. 
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So our jet hadrons finally carry "some" transverse mo- 
nientuixL from fluid partons, but only a small fraction. 
But this is enough to create for example anisotropics in 
dihadron correlations. It will also affect strongly baryon 
to meson rations, as we are going to discuss in a separate 
publication. 



XI. SUMMARY 

We presented a theoretical scheme which accounts for 
bulk matter, jets, and the interaction between the two. 
The criterion for bulk-jet separation is based on parton 
energy loss. But in addition to the latter mechanism, 
there are very important new phenomena which have 
not been discussed so far: The interaction between jet 
hadrons and soft ones (from fluid freeze-out), and the 
interaction between the fluid and jet hadrons at the mo- 
ment of the creation of the latter ones. Particle produc- 
tion between zero and (at least) 20 GeV/c is affected. 
We understand quantitatively azimuthal anisotropics in 



single particle production and dihadron (long range) cor- 
relations at large values of pt . 
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